United States marijuana legalization and opioid mortality trends before and during the first year of the COVID-19 pandemic.

BACKGROUND: To determine if marijuana legalization was associated with reduced opioid mortality. STUDY DESIGN: The United States (US) opioid mortality trend during the 2010-2019 decade was compared in states and District of Columbia (jurisdictions) that had implemented marijuana legalization with states that had not. Acceleration of opioid mortality during 2020, the first year of the coronavirus disease 2019 (COVID-19) pandemic, was also compared in recreational and medicinal-only legalizing jurisdictions. METHODS: Joinpoint methodology was applied to the Centers for Disease Control and Prevention WONDER data. Trends in legalizing jurisdictions were cumulative aggregates. RESULTS: The overall opioid and fentanyl death rates and the percentage of opioid deaths due to fentanyl increased more during 2010-2019 in jurisdictions that legalized marijuana than in those that did not (pairwise comparison p = 0.007, 0.05, and 0.006, respectively). By 2019, the all-opioid and fentanyl death rates were 44 and 50 percent greater in the legalizing than in the nonlegalizing jurisdictions, respectively. When the COVID-19 pandemic hit in 2020, jurisdictions that implemented recreational marijuana legalization before 2019 had significantly greater increases in both overall opioid and fentanyl death rates than jurisdictions with medicinal-only legalization. For all-opioids, the mean (95 percent confidence interval) 2019-to-2020 increases were 46.5 percent (36.6, 56.3 percent) and 29.1 percent (20.2, 37.9 percent), respectively (p = 0.02). For fentanyl, they were 115.6 percent (80.2, 151.6 percent) and 55.4 percent (31.6, 79.2 percent), respectively (p = 0.01). CONCLUSIONS: During the past decade, marijuana legalization in the US was associated at the jurisdiction level with a greater acceleration in opioid death rate. An even greater increase in opioid mortality occurred in recreational-legalizing jurisdictions with the onset of the COVID-19 pandemic. Marijuana legalization is correlated with worsening of the US opioid epidemic.

Pelagic Sargassum in the Gulf of Mexico driven by ocean currents and eddies.

Pelagic Sargassum in the Gulf of Mexico (GoM) plays an important role in ocean biology and ecology, yet our knowledge of its origins and transport pathways is limited. Here, using satellite observations of Sargassum areal density and ocean surface currents between 2000 and 2023, we show that large amounts of Sargassum in the GoM can either originate from the northwestern GoM or be a result of physical transport from the northwestern Caribbean Sea, both with specific transport pathways. Sargassum of the northwestern GoM can be transported to the eastern GoM by ocean currents and eddies, eventually entering the Sargasso Sea. Sargassum entering the GoM from the northwestern Caribbean Sea can be transported in three different directions, with the northward and eastward transports governed by the Loop Current System (LCS) and westward transport driven by the westward extension of the LCS, the propagation or relaying of ocean eddies, the wind-driven westward currents on the Campeche Bank with or without eddies, and the westward currents with/without currents associated with eddies in the northern/central GoM. Overall, the spatial distribution patterns of pelagic Sargassum in the GoM are strongly influenced by the LCS and relevant eddies.

Naturally acquired antibodies against Plasmodium vivax pre-erythrocytic stage vaccine antigens inhibit sporozoite invasion of human hepatocytes in vitro.

In Plasmodium vivax, the most studied vaccine antigens are aimed at blocking merozoite invasion of erythrocytes and disease development. Very few studies have evaluated pre-erythrocytic (PE) stage antigens. The P. vivax circumsporozoite protein (CSP), is considered the leading PE vaccine candidate, but immunity to CSP is short-lived and variant specific. Thus, there is a need to identify other potential candidates to partner with CSP in a multivalent vaccine to protect against infection and disease. We hypothesize that sporozoite antigens important for host cell infection are considered potential targets. In this study, we evaluated the magnitude and quality of naturally acquired antibody responses to four P. vivax PE antigens: sporozoite surface protein 3 (SSP3), sporozoite protein essential for traversal 1 (SPECT1), cell traversal protein of ookinetes and sporozoites (CelTOS) and CSP in plasma of P. vivax infected patients from Thailand. Naturally acquired antibodies to these antigens were prevalent in the study subjects, but with significant differences in magnitude of IgG antibody responses. About 80% of study participants had antibodies to all four antigens and only 2% did not have antibodies to any of the antigens. Most importantly, these antibodies inhibited sporozoite infection of hepatocytes in vitro. Significant variations in magnitude of antigen-specific inhibitory antibody responses were observed with individual samples. The highest inhibitory responses were observed with anti-CelTOS antibodies, followed by anti-SPECT1, SSP3 and CSP antibodies respectively. These data highlight the vaccine potential of these antigens in protecting against hepatocyte infection and the need for a multi-valent pre-erythrocytic vaccine to prevent liver stage development of P. vivax sporozoites.

Interpretable Performance Models for Energetic Materials using Parsimonious Neural Networks.

Predictive models for the performance of explosives and propellants are important for their design, optimization, and safety. Thermochemical codes can predict some of these properties from fundamental quantities such as density and formation energies that can be obtained from first principles. Models that are simpler to evaluate are desirable for efficient, rapid screening of material screening. In addition, interpretable models can provide insight into the physics and chemistry of these materials that could be useful to direct new synthesis. Current state-of-the-art performance models are based on either the parametrization of physics-based expressions or data-driven approaches with minimal interpretability. We use parsimonious neural networks (PNNs) to discover interpretable models for the specific impulse of propellants and detonation velocity and pressure for explosives using data collected from the open literature. A combination of evolutionary optimization with custom neural networks explores and trains models with objective functions that balance accuracy and complexity. For all three properties of interest, we find interpretable models that are Pareto optimal in the accuracy and simplicity space.

Remodelling of Skeletal Muscle Myosin Metabolic States in Hibernating Mammals.

Hibernation is a period of metabolic suppression utilized by many small and large mammal species to survive during winter periods. As the underlying cellular and molecular mechanisms remain incompletely understood, our study aimed to determine whether skeletal muscle myosin and its metabolic efficiency undergo alterations during hibernation to optimize energy utilization. We isolated muscle fibers from small hibernators, Ictidomys tridecemlineatus and Eliomys quercinus and larger hibernators, Ursus arctos and Ursus americanus. We then conducted loaded Mant-ATP chase experiments alongside X-ray diffraction to measure resting myosin dynamics and its ATP demand. In parallel, we performed multiple proteomics analyses. Our results showed a preservation of myosin structure in U. arctos and U. americanus during hibernation, whilst in I. tridecemlineatus and E. quercinus, changes in myosin metabolic states during torpor unexpectedly led to higher levels in energy expenditure of type II, fast-twitch muscle fibers at ambient lab temperatures (20°C). Upon repeating loaded Mant-ATP chase experiments at 8°C (near the body temperature of torpid animals), we found that myosin ATP consumption in type II muscle fibers was reduced by 77-107% during torpor compared to active periods. Additionally, we observed Myh2 hyper-phosphorylation during torpor in I. tridecemilineatus, which was predicted to stabilize the myosin molecule. This may act as a potential molecular mechanism mitigating myosin-associated increases in skeletal muscle energy expenditure during periods of torpor in response to cold exposure. Altogether, we demonstrate that resting myosin is altered in hibernating mammals, contributing to significant changes to the ATP consumption of skeletal muscle. Additionally, we observe that it is further altered in response to cold exposure and highlight myosin as a potentially contributor to skeletal muscle non-shivering thermogenesis.

Climate change is altering the physiology and phenology of an arctic hibernator.

Climate warming is rapid in the Arctic, yet impacts to biological systems are unclear because few long-term studies linking biophysiological processes with environmental conditions exist for this data-poor region. In our study spanning 25 years in the Alaskan Arctic, we demonstrate that climate change is affecting the timing of freeze-thaw cycles in the active layer of permafrost soils and altering the physiology of arctic ground squirrels (Urocitellus parryii). Soil freeze has been delayed and, in response, arctic ground squirrels have delayed when they up-regulate heat production during torpor to prevent freezing. Further, the termination of hibernation in spring has advanced 4 days per decade in females but not males. Continued warming and phenological shifts will alter hibernation energetics, change the seasonal availability of this important prey species, and potentially disrupt intraspecific interactions.

Analysis of heparin-induced thrombocytopenia diagnostic and management strategies in individuals with inconclusive antibody optical densities.

Heparin-induced thrombocytopenia (HIT) is an uncommon but serious complication of exposure to heparin. Antibody optical densities (ODs) used to diagnose HIT exceeding 2 are highly suggestive of disease, whereas ODs less than 0.5 often 'rule out' HIT. Variation in the clinical care of patients with inconclusive ODs between 0.5 and 2 is likely. This single-centre, retrospective analysis evaluates the diagnosis, management and outcomes of those with antibody ODs between 0.5 and 2. We queried our institution's Healthcare Enterprise Repository for Ontological Narration (HERON) database to identify individuals with antibody ODs between 0.5 and 2. Chart review was completed to calculate 4T scores, corroborate diagnosis codes with documented information in our electronic health record (EHR) and evaluate the diagnosis, management and outcomes of these individuals. These data were evaluated using descriptive and univariate statistics. Among individuals evaluated for HIT between November 2007 and July 2020, we identified 302 individuals with ODs between 0.5 and 2. Serotonin release assays (SRAs) were assessed in 55% (165/302) and were positive in 12% (20/165). In those with available data, 96% with low 4T scores had negative SRAs and 4% had positive SRAs. As 4T scores and antibody ODs proportionally increased, SRA positivity also increased. Clinical management varied widely; however, 4T scoring remains a valuable assessment in this cohort. In those with HIT antibody ODs between 0.5 and 2, true positives were uncommon, and their clinical management varied widely. Fortunately, 4T scoring is a useful prognostic tool that improves the diagnosis and management among those with inconclusive HIT.

Risk-Based Monitoring in Clinical Trials: 2021 Update.

Clinical trial quality depends on ensuring participant safety and data integrity, which require careful management throughout the trial lifecycle, from protocol development to final data analysis and submission. Recent developments-including new regulatory requirements, emerging technologies, and trial decentralization-have increased adoption of risk-based monitoring (RBM) and its parent framework, risk-based quality management (RBQM) in clinical trials. The Association of Clinical Research Organizations (ACRO), recognizing the growing importance of these approaches, initiated an ongoing RBM/RBQM landscape survey project in 2019 to track adoption of the eight functional components of RBQM. Here we present results from the third annual survey, which included data from 4889 clinical trials ongoing in 2021. At least one RBQM component was implemented in 88% of trials in the 2021 survey, compared with 77% in 2020 and 53% in 2019. The most frequently implemented components in 2021 were initial and ongoing risk assessments (80 and 78% of trials, respectively). Only 7% of RBQM trials were Phase IV, while the proportions of Phase I-III trials ranged 27-36%. Small trials (< 300 participants) accounted for 60% of those implementing RBQM. The therapeutic areas with the largest number of RBQM trials were oncology (38%), neurology (10%), and infectious diseases (9%). The 2021 survey confirmed a pattern of increasing RBM/RBQM adoption seen in earlier surveys, with risk assessments, which have broad regulatory support, driving RBQM growth; however, one area requiring further development is implementation of centralized monitoring combined with reductions in source data verification (SDV) and source data review (SDR).

A Repeat Pattern of Founder Events for SARS-CoV-2 Variants in Alaska.

Alaska is a unique US state because of its large size, geographically disparate population density, and physical distance from the contiguous United States. Here, we describe a pattern of SARS-CoV-2 variant emergence across Alaska reflective of these differences. Using genomic data, we found that in Alaska, the Omicron sublineage BA.2.3 overtook BA.1.1 by the week of 27 February 2022, reaching 48.5% of sequenced cases. On the contrary, in the contiguous United States, BA.1.1 dominated cases for longer, eventually being displaced by BA.2 sublineages other than BA.2.3. BA.2.3 only reached a prevalence of 10.9% in the contiguous United States. Using phylogenetics, we found evidence of potential origins of the two major clades of BA.2.3 in Alaska and with logistic regression estimated how it emerged and spread throughout the state. The combined evidence is suggestive of founder events in Alaska and is reflective of how Alaska's unique dynamics influence the emergence of SARS-CoV-2 variants.

Pattern of SARS-CoV-2 variant B.1.1.519 emergence in Alaska.

Alaska has the lowest population density in the United States (US) with a mix of urban centers and isolated rural communities. Alaska's distinct population dynamics compared to the contiguous US may have contributed to unique patterns of SARS-CoV-2 variants observed in early 2021. Here we examined 2323 SARS-CoV-2 genomes from Alaska and 278,635 from the contiguous US collected from December 2020 through June 2021 because of the notable emergence and spread of lineage B.1.1.519 in Alaska. We found that B.1.1.519 was consistently detected from late January through June of 2021 in Alaska with a peak prevalence in April of 77.9% unlike the rest of the US at 4.6%. The earlier emergence of B.1.1.519 coincided with a later peak of Alpha (B.1.1.7) compared to the contiguous US. We also observed differences in variant composition over time between the two most populated regions of Alaska and a modest increase in COVID-19 cases during the peak incidence of B.1.1.519. However, it is difficult to disentangle how social dynamics conflated changes in COVID-19 during this time. We suggest that the viral characteristics, such as amino acid substitutions in the spike protein, likely contributed to the unique spread of B.1.1.519 in Alaska.

Building Chemical Property Models for Energetic Materials from Small Datasets Using a Transfer Learning Approach.

For many experimentally measured chemical properties that cannot be directly computed from first-principles, the existing physics-based models do not extrapolate well to out-of-sample molecules, and experimental datasets themselves are too small for traditional machine learning (ML) approaches. To overcome these limitations, we apply a transfer learning approach, whereby we simultaneously train a multi-target regression model on a small number of molecules with experimentally measured values and a large number of molecules with related computed properties. We demonstrate this methodology on predicting the experimentally measured impact sensitivity of energetic crystals, finding that both characteristics of the computed dataset and model architecture are important to prediction accuracy of the small experimental dataset. Our directed-message passing neural network (D-MPNN) ML model using transfer learning outperforms direct-ML and physics-based models on a diverse test set, and the new methods described here are widely applicable to modeling many other structure-property relationships.

Effects of Spring Warming on Seasonal Neuroendocrinology and Activation of the Reproductive Axis in Hibernating Arctic Ground Squirrels.

Many animals adjust the timing of seasonal events, such as reproduction, molt, migration, and hibernation, in response to interannual variation and directional climate-driven changes in temperature. However, the mechanisms by which temperature influences seasonal timing are relatively under-explored. Seasonal timing involves retrograde signaling in which thyrotropin (TSH) in the pars tuberalis (PT) alters expression of thyroid hormone (TH) deiodinases (Dio2/Dio3) in tanycyte cells lining the third ventricle of the hypothalamus. This, in turn, affects the availability of triiodothyronine (T3) within the mediobasal hypothalamus-increased hypothalamic T3 restores a summer phenotype and activates the reproductive axis in long-day breeders. Recently, we showed that retrograde TH signaling is activated during late hibernation in arctic ground squirrels (Urocitellus parryii) held in constant darkness and constant ambient temperature. Sensitivity of seasonal pathways to nonphotic cues, such as temperature, is likely particularly important to hibernating species that are sequestered in hibernacula during spring. To address this issue, we exposed captive arctic ground squirrels of both sexes to an ecologically relevant increase in ambient temperature (from -6 to -1°C) late in hibernation and examined the effects of warming on the seasonal retrograde TSH/Dio/T3 signaling pathway, as well as downstream elements of the reproductive axis. We found that warmed males tended to have higher PT TSHß expression and significantly heavier testis mass whereas the TSH/Dio/T3 signaling pathway was unaffected by warming in females, although warmed females exhibited a slight decrease in ovarian mass. Our findings suggest that temperature could have different effects on gonadal growth in male and female arctic ground squirrels, which could lead to mismatched timing in response to rapid climate change.

A repeat pattern of founder events for SARS-CoV-2 variants in Alaska.

Alaska is a unique US state because of its large size, geographically disparate population density, and physical distance from the contiguous United States. Here, we describe a pattern of SARS-CoV-2 variant emergence across Alaska reflective of these differences. Using genomic data, we found that in Alaska the Omicron sublineage BA.2.3 overtook BA.1.1 by the week of 2022-02-27, reaching 48.5% of sequenced cases. On the contrary in the contiguous United States, BA.1.1 dominated cases for longer, eventually being displaced by BA.2 sublineages other than BA.2.3. BA.2.3 only reached a prevalence of 10.9% in the contiguous United States. Using phylogenetics, we found evidence of potential origins of the two major clades of BA.2.3 in Alaska and with logistic regression estimated how it emerged and spread throughout the state. The combined evidence is suggestive of founder events in Alaska and is reflective of how Alaska’s unique dynamics influence the emergence of SARS-CoV-2 variants.

Insulin-like growth factor-I biocompartmentalization across blood, interstitial fluid and muscle, before and after 3 months of chronic resistance exercise.

This investigation examined the influence of 12-week ballistic resistance training programs on the IGF-I system in circulation, interstitial fluid, and skeletal muscle, at rest and in response to acute exercise. Seventeen college-aged subjects (11 women/6 men; 21.7 ± 3.7 yr) completed an acute ballistic exercise bout before and after the training program. Blood samples were collected pre-, mid-, and postexercise and analyzed for serum total IGF-I, free IGF-I, and IGF binding proteins (IGFBPs) 1-4. Dialysate and interstitial free IGF-I were analyzed in vastus lateralis (VL) interstitial fluid collected pre- and postexercise via microdialysis. Pre- and postexercise VL muscle biopsies were analyzed for IGF-I protein expression, IGF-I receptor phosphorylation (p-IGF-IR), and AKT phosphorylation (p-AKT). Following training, basal serum IGF-I, free IGF-I, IGFBP-2, and IGFBP-3 decreased whereas IGFBP-1 and IGFBP-4 increased. Training reduced basal dialysate and interstitial free IGF-I but had no effect on basal skeletal muscle IGF-I, p-IGF-IR, or p-AKT. Acute exercise elicited transient changes in IGF-I system concentrations and downstream anabolic signaling both pre- and posttraining; training did not affect this acute exercise response. Posttraining, acute exercise-induced changes in dialysate/interstitial free IGF-I were strongly correlated with the changes in intramuscular IGF-I expression, p-IGF-IR, and p-AKT. The divergent influence of resistance training on circulating/interstitial and skeletal muscle IGF-I demonstrates the importance of concurrent, multiple biocompartment analysis when examining the IGF-I system. As training elicited muscle hypertrophy, these findings indicate that IGF-I's anabolic effects on skeletal muscle are mediated by local, rather than systemic mechanisms.NEW & NOTEWORTHY In the first investigation to assess resistance training's effects on the IGF-I system in serum, interstitial fluid, and skeletal muscle, training decreased basal circulating and interstitial IGF-I but did not alter basal intramuscular IGF-I protein activity. Posttraining, acute exercise-induced interstitial IGF-I increases were strongly correlated with intramuscular IGF-I expression and signaling. These findings highlight the importance of multibiocompartment measurement when analyzing IGF-I and suggest that IGF-I's role in hypertrophic adaptations is locally mediated.

Scale dependence of coral reef oases and their environmental correlates.

Identifying relatively intact areas within ecosystems and determining the conditions favoring their existence is necessary for effective management in the context of widespread environmental degradation. In this study, we used 3766 surveys of randomly selected sites in the United States and U.S. Territories to identify the correlates of sites categorized as "oases" (defined as sites with relatively high total coral cover). We used occupancy models to evaluate the influence of 10 environmental predictors on the probability that an area (21.2-km2 cell) would harbor coral oases defined at four spatial extents: cross-basin, basin, region, and subregion. Across all four spatial extents, oases were more likely to occur in habitats with high light attenuation. The influence of the other environmental predictors on the probability of oasis occurrence were less consistent and varied with the scale of observation. Oases were most likely in areas of low human population density, but this effect was evident only at the cross-basin and subregional extents. At the regional and subregional extents oases were more likely where sea-surface temperature was more variable, whereas at the larger spatial extents the opposite was true. By identifying the correlates of oasis occurrence, the model can inform the prioritization of reef areas for management. Areas with biophysical conditions that confer corals with physiological resilience, as well as limited human impacts, likely support coral reef oases across spatial extents. Our approach is widely applicable to the development of conservation strategies to protect biodiversity and ecosystems in an era of magnified human disturbance.

Hypothalamic remodeling of thyroid hormone signaling during hibernation in the arctic ground squirrel.

Hibernation involves prolonged intervals of profound metabolic suppression periodically interrupted by brief arousals to euthermy, the function of which is unknown. Annual cycles in mammals are timed by a photoperiodically-regulated thyroid-hormone-dependent mechanism in hypothalamic tanycytes, driven by thyrotropin (TSH) in the pars tuberalis (PT), which regulates local TH-converting deiodinases and triggers remodeling of neuroendocrine pathways. We demonstrate that over the course of hibernation in continuous darkness, arctic ground squirrels (Urocitellus parryii) up-regulate the retrograde TSH/Deiodinase/TH pathway, remodel hypothalamic tanycytes, and activate the reproductive axis. Forcing the premature termination of hibernation by warming animals induced hypothalamic deiodinase expression and the accumulation of secretory granules in PT thyrotrophs and pituitary gonadotrophs, but did not further activate the reproductive axis. We suggest that periodic arousals may allow for the transient activation of hypothalamic thyroid hormone signaling, cellular remodeling, and re-programming of brain circuits in preparation for the short Arctic summer.

United States marijuana legalization and opioid mortality epidemic during 2010-2020 and pandemic implications.

BACKGROUND: The hypothesis that marijuana availability reduces opioid mortality merits more complete testing, especially in a country with the world's highest opioid death rate and 2nd highest cannabis-use-disorder prevalence. METHODS: The United States opioid mortality rate was compared in states and District of Columbia that had implemented marijuana legalization with states that had not, by applying joinpoint methodology to Centers for Disease Control and Prevention data. Variables included race/ethnicity and fentanyl-type opioids (fentanyls). RESULTS: After the same rates during 2010-2012, the opioid mortality rate increased more rapidly in marijuana-legalizing than non-legalizing jurisdictions (2010-2020 annual pairwise comparison p = 0.003 for all opioids and p = 0.0004 for fentanyls). During the past decade, all four major race/ethnicities in the U.S. had evidence for a statistically-significant greater increase in opioid mortality rates in legalizing than non-legalizing jurisdictions. Among legalizing jurisdictions, the greatest mortality rate increase for all opioids was in non-Hispanic blacks (27%/year, p = 0.0001) and for fentanyls in Hispanics (45%/year, p = 0.0000008). The greatest annual opioid mortality increase occurred in 2020, the first year of the COVID-19 pandemic, with non-Hispanic blacks having the greatest increase in legalizing vs. non-legalizing opioid-death-rate difference, from 32% higher in legalizing jurisdictions in 2019 to more than double in 2020. CONCLUSIONS: Instead of supporting the marijuana protection hypothesis, ecologic associations at the national level suggest that marijuana legalization has contributed to the U.S.'s opioid epidemic in all major races/ethnicities, and especially in blacks. If so, the increased use of marijuana during the 2020-2022 pandemic may thereby worsen the country's opioid crisis.

Risk-Based Monitoring in Clinical Trials: Increased Adoption Throughout 2020.

With the emergence of new technologies for data collection, the continued impact of the COVID-19 pandemic, and the increasing number of partially or fully decentralized clinical trials (DCTs), the importance of risk-based monitoring (RBM) and the larger risk-based quality management (RBQM) framework in clinical trial management is increasing. RBM and RBQM focus on the detection of events or trends that impact trial quality in terms of participant safety and data integrity. In 2019, the Association of Clinical Research Organizations (ACRO) began a landscape survey of RBM/RBQM implementation in ongoing clinical trials. Initial results of this survey, representing full-year data for 2019, were reported previously. Here, we present full-year landscape data for 2020 drawn from 5,987 clinical trials ongoing at the end of 2020, including 908 new studies started that year. Of these trials, 77% implemented at least one RBM/RBQM component, an increase from 47% for studies ongoing at the end of 2019. We also observed increased implementation for three of the five RBM components included in the survey. Centralized monitoring decreased nominally in 2020 compared with 2019. Although the percentages of 2020 trials incorporating reduced source data verification (SDV) and reduced source data review (SDR) increased from 2019 to 2020, these numbers are still low considering the large percentage of trials implementing at least one RBQM component. In the current clinical trial landscape, as more DCTs are launched and new data collection technologies are implemented, there remains a pressing need for greater use of centralized monitoring coupled with reductions in SDR/SDV and, ultimately, greater adoption of RBM and RBQM.

Transcriptional changes and preservation of bone mass in hibernating black bears.

Physical inactivity leads to losses of bone mass and strength in most mammalian species. In contrast, hibernating bears show no bone loss over the prolonged periods (4-6 months) of immobility during winter, which suggests that they have adaptive mechanisms to preserve bone mass. To identify transcriptional changes that underlie molecular mechanisms preventing disuse osteoporosis, we conducted a large-scale gene expression screening in the trabecular bone and bone marrow, comparing hibernating and summer active bears through sequencing of the transcriptome. Gene set enrichment analysis showed a coordinated down-regulation of genes involved in bone resorption, osteoclast differentiation and signaling, and apoptosis during hibernation. These findings are consistent with previous histological findings and likely contribute to the preservation of bone during the immobility of hibernation. In contrast, no significant enrichment indicating directional changes in gene expression was detected in the gene sets of bone formation and osteoblast signaling in hibernating bears. Additionally, we revealed significant and coordinated transcriptional induction of gene sets involved in aerobic energy production including fatty acid beta oxidation, tricarboxylic acid cycle, oxidative phosphorylation, and mitochondrial metabolism. Mitochondrial oxidation was likely up-regulated by transcriptionally induced AMPK/PGC1α pathway, an upstream stimulator of mitochondrial function.